Designing Topological Bands in Reciprocal Space 
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Motivated by new capabilities to realise artificial gauge fields in ultracold atomic systems, and 
by their potential to access correlated topological phases in lattice systems, we present a new strat- 
egy for designing topologically non-trivial band structures. Our approach is simple and direct: it 
amounts to considering tight-binding models directly in reciprocal space. These models naturally 
cause atoms to experience highly uniform magnetic flux density and lead to topological bands with 
very narrow dispersion, without fine-tuning of parameters. Further, our construction immediately 
yields instances of optical Chern lattices, as well as band structures of higher Chern number, \C\ > 1. 



The ability to create new environments for custom- 
tailored degrees of freedom is central to advancing our 
understanding of correlated electron physics. In the same 
way that modulation-doped semiconductor heterostruc- 
tures gave rise to the field of fractional quantum Hall 
physics, the creation of artifical gauge fieldsp] holds 
the promise of providing an entirely new perspective 
on strong-correlation physics in topological bands. The 
question that naturally arises now is not only how to gen- 
erate such artificial gauge fields, but also what settings 
will be most promising for realising new physical phenom- 
ena. There is intense current interest in the properties 
of such systems in condensed matter physics, most re- 
cently through the proposed existence of fractional Chern 
insulators [2HS]. These are lattice systems in which frac- 
tional quantum Hall physics occurs in partially filled non- 
dispersive topological "Chern" bands O [7]. These works 
are all based on tight-binding models, in which band 
topology and dispersion are controlled by the tunneling 
matrix elements. The first steps have been taken in im- 
plementing tight-binding models of this type for atoms 
in optical lattices [8j, with tunable tunneling phases [9] at 
least for nearest-neighbour sites [TUHH]. 

Recently it has been demonstrated that topological 
bands can be formed for ultracold atoms in another, very 
direct, wav |13Ifl5] . Coherent Raman coupling of several 
(N > 1) internal atomic statespQ, as used to generate 
artificial gauge fields[TH], allows new forms of optical lat- 
tice. In particular, in "optical flux lattices" atoms expe- 
rience non-zero average magnetic flux density, resulting 
in low-energy topological bands that are analogous to 
Landau levels [TBI ITS] . These new forms of optical lattice 
are readily implemented in experiment, so it is highly 
desirable to be able to design their properties. However, 
they can operate far from the tight-binding limit, so the 
existing solid-state models cannot be directly applied. 

Here, we present a number of insights into the nature 
of optical flux lattices and develop theoretical tools that 
allow a targeted design of topological bands in optical lat- 
tices. Notably, we show that optical flux lattices are use- 
fully viewed as tight-binding models in reciprocal space. 
This allows a design of band structures by transposing 



established results on Chern bands in real space into our 
new setting. We use our insights to construct three new 
forms of optical lattice. First, we show how to generate 
optical flux lattices, involving more than N — 2 inter- 
nal states, for which the magnetic flux density is highly 
uniform in space. The resulting energy bands have very 
narrow bandwidth, and recover continuum Landau levels 
for large N. Second, we show how to construct "optical 
Chern lattices", in which the dressed states experience 
zero net magnetic flux, but for which the band has non- 
zero Chern number. Finally, we show how one can de- 
sign optical lattices in which the lowest energy band has 
Chern number larger than one. 

Our results follow from examining the limits of 
deep/shallow lattices, when the lattice depth V is 
large/small compared to the typical kinetic energy Er = 
h 2 K 2 /(2M) (all symbols are defined below). After setting 
up the problem, we discuss these two limits in turn. We 
comment on possible implementations on the way. 

The Hamiltonian for an optically dressed atom with N 
internal states isp] 

H=^i N + V(r) (1) 

where Ijv is the N x N identity matrix, and V(r) an 
N x N matrix describing the coherent optical coupling 
in the rotating wave approximation. The nature of the 
internal states is unimportant for our purposes: these can 
be electronic excited states, hyperfine states, or states of 
orbital motion. 

In all cases, the coupling V(r) involves absorp- 
tion/emission of photons from laser beams, so it natu- 
rally appears as couplings Vff" k = (a' , fe'|y|a, k), where 
\a,k) is the momentum k eigenstate for component a. 
We focus on periodic lattices, such that the set of mo- 
mentum transfers {k,} for which V£ a is non-zero forms 
a regular lattice. The momentum k of any given com- 
ponent a is then only conserved up to the addition of 
reciprocal lattice vectors, the basis vectors of which we 
denote Gi, with i = 1 . . . d where d is the dimensionality 
of the lattice. (We focus on lattices in d = 2, but the 
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results can be extended to d — 3, with Berry curvature 
and magnetic flux density becoming pseudovectors.) 

By Bloch's theorem, the energy eigenstates of ([I]) can 
be assigned crystal momentum q and band index n, and 
be decomposed as \ip nq ) = I] a>G c"£,|a!, g - g a - G) 
where a runs over all N internal states, G runs over all 
sites of the reciprocal lattice, and the vectors g a account 
for possible momentum offsets involved in the inter-state 
transitions. The energy eigenvalues E nq follow from 
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where e k = h 2 \k\ 2 /(2M). 

Deep lattice/ adiabatic limit: In the limit £r/V — > the 
kinetic energy can be neglected. The effects of the opti- 
cal coupling in Eqn. can then be fully understood in 
terms of the local dressed states, the eigenstates of V(r). 
Under adiabatic motion, the n th dressed state experi- 
ences magnetic flux density [T] n^r) = J2 a Vu™* x 
Vw™, due to the Berry curvature [17] associated with spa- 
tial variations of its wavefunction u a (r). In an optical 
flux lattice, the lowest energy dressed state experiences 
7^ magnetic flux quanta through each unit cell [13]. 

It is instructive to examine this limit -Er/V — > also 
from the point of view of Eqn. Dropping the kinetic 
energy causes this eigenvalue problem to reduce to that 
of a uniform tight-binding model in reciprocal space, de- 
fined by the couplings V^" G _ fl f _ G , between sites at 
positions g a + G. The energy spectrum of this reciprocal- 
space tight-binding model - its "bandstructure" - con- 
sists of N bands, since there are N sites g a associated 
with each reciprocal lattice point G. The eigenstates 
are extended Bloch waves in reciprocal space, and can 
be written c n J G cx e li - G+ ^ r u n a (r). (For E R /V -> the 
energy is independent of q and the states are naturally 
labelled by the band index n = 1 . . . N and a conserved 
"momentum" r.) 

Crucially, we identify this conserved "momentum" of 
the reciprocal-space tight-binding model with the real 
space position r (and thus its "Brillouin zone" with the 
real space unit cell). The state in the n th band sim- 
ply corresponds to the n th local dressed state of V(r). 
The power of this viewpoint emerges upon considering 
the magnetic flux density n^: the Berry curvature asso- 
ciated with the adiabatic motion of any dressed state in 



real space equals the Berry curvature of the associated 
band of the reciprocal-space tight-binding model. In par- 
ticular, the number of magnetic flux quanta through 
the real space unit cell equals the Chern number of that 
band. Hence, the criterion for an optical flux lattice is 
that the lowest energy band of the reciprocal- space tight- 
binding model has non-zero Chern number. Since much 
is known about the Chern bands of various tight-binding 
models, this criterion can be used to design optical flux 
lattices with specific properties. 




FIG. 1: One primitive unit cell of the reciprocal-space tight- 
binding model corresponding to the triangular lattice optical 
flux lattice with TV internal states J3j, showing the optical 
couplings — Ve 1 ^ . The clockwise sum of the phases <f> around 
any one of the triangular plaquettes is n/N modulo 2tt. 



One highly desirable feature is to generate a mag- 
netic flux density that is uniform. The resulting energy 
bands will then closely reproduce Landau levels: topo- 
logical bands with exact degeneracy, highly susceptible 
to strong-correlation physics. For all previously proposed 
optical flux lattices [TTJ [T31 [T5] the magnetic flux density 
is very non-uniform, vanishing at a set of points. This 
non- uniformity is a mathematical necessity for N = 2|13j 
and for general N when V consists of the generators of 
SU(2), as in Ref. The above considerations show how 
to overcome this limitation: one should form a reciprocal- 
space tight-binding model for which the lowest energy 
band has uniform Berry curvature. 

An optical flux lattice that achieves this goal is illus- 
trated in Fig. [TJ The N internal states are arranged on 
the sites of a triangular lattice, marked by the state la- 
bel a. The links between lattice sites indicate optical 
couplings, the displacement in reciprocal space the mo- 
mentum transfer, with the label marking amplitude and 
phase. Explicitly, Fig. [l] encodes the coupling matrix 
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FIG. 2: Spatial dependence of (a) the energy and (b) the 
magnetic flux density of the lowest energy dressed state of 
the triangular optical flux lattice (|3| for N = 4, within a 
unit cell containing N# = 1 flux quantum. Both are highly 
uniform, becoming increasingly so as AT — > oo. 



where Aj = exp(— ir ■ Kj), Ki = (1, 0)k , «2 = ( \, 2 ) K 
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and = Ki — «2- We choose the primitive unit cell of 
the reciprocal lattice to have basis vectors G\ = Nk% 
and Gi = K3; the real-space unit cell then has lattice 
vectors o x = ^(V5/2, 1/2) and a 2 = ^(0, 1). 

The properties of the resulting dressed states follow 
from the bandstructure of the reciprocal-space tight- 
binding model. The phases on the nearest-neighbour 
couplings mimic the effect of a uniform magnetic field 
in reciprocal space (the total phase acquired on hopping 
around any triangular plaquette being tt/N, modulo 2n), 
so the spectrum follows from Harper's equation for the 
triangular lattice [TS]. For N > 2 the lowest energy band 
has Chern number C = 1, implying an optical flux lattice 
with N$ = 1 flux quantum per unit cell, i.e. mean flux 
density 
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The local magnetic flux density and the dressed state 
energy follow from the Berry curvature and dispersion 
of the lowest energy band of the reciprocal-space tight- 
binding model. Both show spatial modulations, illus- 
trated in Fig. [2] for N — 4. (The spatial dependence 
of the dressed state takes a form akin to a triangular 
multi-component Skyrmion lattice [TfJ] with lattice con- 
stant |ai|.) However, as N — > 00, the reciprocal-space 
tight-binding model recovers the continuum limit, and 
its lowest band becomes a degenerate lowest Landau level 
with uniform Berry curvature. The convergence to uni- 
formity with increasing N is very fast, with corrections 
that fall exponentially auicklv [T%l ITU] . 

The emergence, for large N, of highly uniform mag- 
netic flux density distinguishes this optical flux lattice 
from all previous proposals. Owing to this uniformity, 
the resulting bandstructure for non-zero kinetic energy, 
Er V, closely reproduces the spectrum of a charged 
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FIG. 3: Low energy density of states for: (a) The triangular 
optical flux lattice (J3J) with N = 4 internal states and V = Eb_; 
the spectrum closely reproduces the Landau level spectrum, 
all bands shown having Chern number C = 1. (b) A related 
model (still with = 4 and V = £?rt) in which the phases on 
the optical couplings are modified as described in the text to 
form a lowest energy band with Chern number C = 2. 



particle in a uniform magnetic field: degenerate Lan- 
dau levels, spaced by the cyclotron energy Huj c 
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Sr. The fast convergence means that this qualita- 



tive structure is already apparent for small N. This is 
clearly seen in Fig. (3^a) for N = 4, for which the lowest 
energy band has a bandwidth that is about 200 times 
smaller than the gap to the next band. (For N = 3 the 
same ratio is about 40.) 

One can envisage a variety of experimental implemen- 
tations for lattices of this type, or the many variants 
that achieve the same end. The case N = 2 is the 
triangular optical flux lattice described in Ref. [T3j For 
N = 3 one can use the three components of the F = 1 
hyperfine manifold of 87 RbJ5U]. A simple variant is ob- 
tained by omitting the state-dependent coupling along 
K3 (the geometry can then be made square by rotating 
K2 — > (0, 1)k). This lattice requires only cyclic coupling 
over the N — > N ± 1 states. An implementation of cou- 
plings of this kind has been proposed for N = 4 using 
hyperfine levels of 87 Rb. 

Weak-lattice limit: Further useful results can be ob- 
tained by considering the eigenvalue problem ^ in the 
"nearly- free-electron" limit V <C Here, the only in- 
fluence of the couplings V£ a is at values of crystal mo- 
mentum q close to the lines where two free-particle states 
are degenerate, where they open up band-gaps. Close to 
those points q* where two (or more) lines cross, the bands 
can also acquire non-zero Berry curvature. Therefore, 
one can determine the Chern numbers of the bands in 
the weak-lattice limit by computing the Berry curvature 
in the vicinity of these high symmetry points q* . 

To illustrate the construction, consider the triangu- 
lar optical flux lattice, Fig. [T] in this limit. The high 
symmetry points, q* , are located at the centres of the 
triangular plaquettes, at which the kinetic energies of 
three states (labelled a, b, c) are equal. (These points are 
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marked by *'s in Fig. [I]) This degeneracy is split by the 
optical couplings around this plaquette, which we denote 
V ab = -Ve**">, V bc = -Ve^«, V ca = -Ve**~. lute- 
grating the Berry curvature for the lowest energy state 
over the vicinity of q* shows that the total contribution is 
0tot = mod (ir + <f) ab + <f) bc + (f> cai 2ir)-ir. (For |^ to t| = w 
the lowest two bands cross at a Dirac point, so the Berry 
curvature of the lowest energy band is not defined.) For 
the triangular lattice model ^ the phases are such that 
0tot — "x/N for each of the 2N plaquettes. Thus we im- 
mediately establish that the Chern number of the lowest 
energy band is ^ x (2N) x (ir/N) = 1. 

Our approach allows us to design other forms of optical 
lattice, with topological bands that are not directly re- 
lated to Landau levels. First we show that one can form 
optical lattices which have bands with non-zero Chern 
number, but for which the adiabatic dressed state in real 
space feels vanishing net flux. In the solid state literature, 
these are referred to as Chern bands. All previous exam- 
ples involve tight-binding models in real space, following 
the seminal work of Haldane [7 . The analogous "opti- 
cal Chern lattices" that we describe here have the new 
feature that the magnetic flux density is defined as a con- 
tinuous function of real-space position. A simple example 
of such an optical Chern lattice is obtained by choosing 
the triangular optical flux lattice, Eqn. ([3]), with N = 2 
and introducing a state-dependent potential SVa z . For 
SV > 2V the lowest energy dressed state has (a z ) < 0. 
Hence, the Bloch vector[13] cannot wrap the sphere - 
the net flux through any unit cell must vanish. However, 
this additional term does not affect the Berry curvatures 
at the (four) symmetry points in the Brillouin zone in 
the weak-lattice limit: the Chern number of the lowest 
energy band remains C = 1. 

Finally, we show how to design an optical lattice for 
which the lowest energy band has Chern number of mag- 
nitude larger than one, \C\ > 1, in the weak-lattice limit. 
One should arrange that the integrated Berry curvatures 
from all symmetry points q* sum to the appropriate value 
2itC. As a concrete example, we adapt the triangular lat- 
tice ([3| to form a lowest energy band with arbitrary C by 
replacing the phases jfr appearing in ^ , in both diagonal 
and off-diagonal entries, by Cjj. The requirement that 
the integrated Berry curvature per plaquette lies in the 
range — it < Cfj < n, leads to \C\ < N. Hence, a lowest 
energy band with C = 2 can be achieved for N > 3 inter- 
nal states. Numerical studies for N = 3,4 confirm that 
the lowest energy band retains C — 2 for V ~ Er, and 
even into the deep-lattice limit. The low energy band- 
structure for N = 4 for V = Er is shown in Fig. [3jb). 
The lowest energy band has Chern number C = 2. The 
band is well separated from higher bands, the gap to the 
next band being about 4.5 times its bandwidth, despite 
the fact that no optimization has been applied. 

The methods we have described open the door for the 
construction of many other forms of optical flux lattice. 



Any tight-binding model for which the lowest energy 
band has non-zero Chern number can be used as the basis 
for an optical flux lattice. Hence uniform magnetic flux 
density in real space could also be achieved for a system 
with fixed N by introducing further-neighbour couplings 
on the reciprocal-space lattice [22]. Such couplings are 
readily implemented as higher momentum transfers; in- 
deed, the optical flux lattice proposed in Ref . [TBI involves 
all the couplings of the Haldane model|7|. We have also 
provided a direct method by which to construct optical 
lattices with low energy Chern bands of any topology in 
the weak-lattice limit. This ability to design topological 
bandstructures, with narrow bandwidth, will allow the 
use of atomic gases to explore novel correlated topologi- 
cal phases of bosons or fermions in lattice systems. 

NRC is grateful to Jean Dalibard for numerous help- 
ful discussions, and RM to Benoit Dougot and Dima 
Kovrizhin for collaboration on Ref. [19]. This work was 
supported by EPSRC Grant EP/F032773/1. We thank 
the Galileo Galilei Institute for hospitality. 



[1] 

[2] 
[3] 

[4: 

[5 
[6 

in 

[8] 
[9 

[10 

[11 

[12 
[13 
[14 

[15 

[16 
[17 

[18 

[19 

[20 
[21 

[22] 



J. Dalibard, F. Gerbier, G. Juzeliunas, and P. Ohberg, 

Rev. Mod. Phys. 83, 1523 (2011). 

R. Roy and S. L. Sondhi, Physics 4, 46 (2011). 

E. Tang, J.-W. Mei, and X.-G. Wen, Phys. Rev. Lett. 
106, 236802 (2011). 

K. Sun, Z. Gu, H. Katsura, and S. Das Sarma, Phys. 
Rev. Lett. 106, 236803 (2011). 

T. Neupert, L. Santos, C. Chamon, and C. Mudry, Phys. 
Rev. Lett. 106, 236804 (2011). 

D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. 
den Nijs, Phys. Rev. Lett. 49, 405 (1982). 

F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
M. Lewenstein et al., Advances in Physics 56, 243 (2007). 
D. Jaksch and P. Zoller, New Journal of Physics 5, 56 
(2003). 

M. Aidelsburger et al, Phys. Rev. Lett. 107, 255301 

(2011) . 

K. Jimenez-Garcfa et al, Phys. Rev. Lett. 108, 225303 

(2012) . 

J. Struck et al, Phys. Rev. Lett. 108, 225304 (2012). 
N. R. Cooper, Phys. Rev. Lett. 106, 175301 (2011). 
B. Beri and N. R. Cooper, Phys. Rev. Lett. 107, 145301 
(2011). 

N. R. Cooper and J. Dalibard, Europhysics Letters 95, 
66004 (2011). 

Y.-J. Lin et al, Nature 462, 628 (2009). 

D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 

1959 (2010). 

F. H. Claro and G. H. Wannier, Phys. Rev. B 19, 6068 
(1979). 

D. L. Kovrizhin, B. Doucot, and R. Moessner, 
larXiv:1207.402T1 

J. Dalibard, private communication. 

D. L. Campbell, G. Juzeliunas, and I. B. Spielman, Phys. 
Rev. A 84, 025602 (2011). 

E. Kapit and E. Mueller, Phys. Rev. Lett. 105, 215303 
(2010). 



